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We study antimony doping effects in the iron-based superconductor CaFe(Sb;cAsi_;c )2 by using the first-principle 
calculation. The calculations reveal that the substitution of the doped antimony atom into As of the chainlike As layers 
is more stable than that in FeAs layers. This prediction can be checked by experiments. Our results suggest that doping 
homologous elements into the chainlike As layers existing only in novel 112 system is responsible for rising up the 
critical temperature. We discuss antimony doping effects on the electronic structure. It is found that the calculated band 
structures with and without the antimony doping are similar to each other within our framework. 


The discovery of the iron-based superconductor 
LaFeAsOi-jcFjc has attracted much attention because of 
the high superconducting transition temperature.^^ Many dif¬ 
ferent families of iron-based materials have been found, such 
as RE¥qA sO (1111 type),^^ AE¥q 2^^2 (122 type),^^ AFeAs 
(111 type),^^ and FeSe (11 type)^^ and so on, which contain 
rare-each (RE), alkali-earth {AE), and alkali (A) elements. 
The maximum superconducting transition temperature is 
55K and all the iron-based superconductors consist of FeAs 
(FeSe) layers and various kinds of blocking layers. 

Recently, the discovery of the novel 112-type iron arsenides 
of Cai_;cLajcFeAs 2 , Cai_;cPfjcFeAs 2 and C 2 ii-xRex¥eAs 2 
(RE = Ce, Nd, Sm, Eu, and Gd) have received considerable 
attention.The crystal structure of the 112-type iron ar¬ 
senides is monoclinic and consists of alternately stacked FeAs 
and arsenic zigzag bond layers. Very recently, it was reported 
that the simultaneous doping of isovalent P or Sb drastically 
improves the superconductivity in the La-doped 112 phase. 
In addition, a large amount of Sb doping leads to a further in¬ 
crease in Tc to 47 Kubo et al. found that the enhancement 

of superconductivity originates from the enlargement in the 
lattice parameter b, which increases the As-Fe-As bond an¬ 
gle to be closer to the ideal tetrahedron value. We note that 
the Tc increase by the antimony substitution does not occur 
in any other iron pnictides. In order to understand the mech¬ 
anism of the enhancement of Tc, we need to know where the 
antimony is substituted, since there are two possible Sb sub¬ 
stitution points; the chainlike As layer and the FeAs layer. 

The electronic structure of CaFeAs 2 was investigated 
by two groups using the first-principle calculations.^’ 
Katayama et al. reported, with the use of the experimental 
lattice parameters, that the Fermi surfaces consist of cylindri¬ 
cal ones and the overall characteristics are similar to those of 
the high-temperature iron-based superconductors LaFeAsO. 
The band calculation by Wu et al. with the optimized lattice 
parameters suggested that the Fermi surfaces are very similar 
to the normal Fermi surfaces in other iron-based superconduc¬ 



Eig. 1. (Color online)2 x 2 x 1 supercell of CaFeAs 2 with Sb substitution 
in (a) the chainlike As layer and (b) the FeAs layer 


tors except for the additional 3D hole pocket and four electron 
cones. 

In this paper, we perform first-principle calculations to in¬ 
vestigate antimony doping effects on CaFe(Sb;cAsi_;c) 2 - Con¬ 
sequently, we find that the antimony tends to be substituted 
into the chainlike As layers, not the FeAs layers, since the 
energy of the Sb substitution into the chainlike As layers is 
0.63eV per 2x2x1 superceh of CaFeAs 2 lower than that 
into the FeAs layer. In these calculations, we optimize the 
lattice parameters and atom positions. In the case of the an¬ 
timony substitution into the chainlike As layer, we observe 
an increase in the lattice parameter b. This result is consistent 
with the experiments. We also show the difference of the band 
dispersions with and without the Sb substitution. 

Here, we briefly show the present calculation techniques. 
The calculation package employed throughout this paper is 
VASP,^^^ which adopts PAW method^^^ and GGA exchange- 
correlation energy. In order to examine where the an¬ 
timony is substituted, we prepare 2x2x1 superceh of 
CagFegSbAsis as shown in Fig. 1. The chemical formula is 
given by CaFe(Sb;cAsi_jc) 2 , in which x = 0.0625. We opti- 
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mize the crystal structure of CaFe(SbjcAsi_;c )2 until the forces 
on all atoms become less than 0.01 eV/A. We relax atom po¬ 
sitions, unit cell shape, and unit cell volume. The parameter 
set in the DFT calculations are given as follows. A:-points are 
taken as 10x11x8, and the self-consistent loops are repeated 
until the energy deviation becomes less than 10“^eV with the 
cut-off energy being 600 eV. 

Now, let us discuss the calculated crystal structure of 
CaFe(Sbo.o 625 Aso. 9375 ) 2 . The obtained lattice parameters are 
summarized in Table I. The parameters in the experimental 
structure for Cai_jcLa;cFeAs 2 was reported in Ref. 5. It should 
be noted that the lattice parameter c in the optimized crystal 
structure is smaller than that obtained by experiments, which 
is usually found in iron-based materials. This might be orig¬ 
inated from the fact that the Van der Waals force can not be 
accurately taken into consideration by the DFT calculation in 
the layer-type materials. We find that the antimony tends to be 
substituted into the chainlike As layer, since the calculated en¬ 
ergy of the antimony substitution into the chainlike As layers 
is 0.63eV per 2x2x1 supercell of CaFeAs 2 lower than that 
into the FeAs layer. Focusing on the crystal structure with the 
antimony substitution into the chainlike As layer, the lattice 
parameters a and h increases with the antimony doping. The 
calculation suggests that the unit volume also increases with 
the antimony doping, since the value of the c axis in the anti¬ 
mony doped system is similar to that in the optimized system. 
The calculations reveal that the substitution of the doped an¬ 
timony atom into the chainlike As layers is more stable than 
that in FeAs layers. This prediction can be checked by exper¬ 
iments. Our results suggest that the chainlike As layers, exist¬ 
ing only in novel 112 system, have crucial role for increasing 
the critical temperature. 

Phenomenologically, it is known that the of the iron- 
based superconductors is optimized around the anion height 
of 1.38 A, or the As-Fe-As bond angle of 109°.^^’^^^ From 
this viewpoint, increase may be expected more when anti¬ 
mony substitution occurs in the chainlike As layer, because in 
this case the height decreases from 1.418 A toward 1.38 A, 
and the angle increases from 108.079° toward 109°. However, 
there are exceptions for these phenomenological rules of op¬ 
timizing Tc, and it is not certain whether the rules apply to the 
present material. Therefore, we need to see how these changes 
in the lattice parameters affect the electronic structure in the 
following. 

We first study the Fermi surfaces for the parent compound 
CaFeAs 2 . We use the structural parameters obtained by the 
experiment in Ref. 5. The calculation results show that there 
are three two-dimensional cylindrical Fermi surfaces and one 
three-dimensional ellipsoidal Fermi surface around F point, 
two cylindrical Fermi surfaces around A = (1/2,1/2,0) 
point, and one very narrow cylindrical Fermi surfaces near 
B = (1/2,0,0) point, as shown in Fig. 2. The Fermi sur¬ 
faces are similar to those in 1111-type iron-based supercon¬ 
ductors, except for the three-dimensional ellipsoidal and very 
narrow cylindrical Fermi surfaces. Note that Wu et al. showed 
the Fermi surfaces with the optimized lattice parameters in 


Table II. Anion height from Fe layer and As-Fe-As bond angle in calcu¬ 
lated crystal structures of CaFe(Sbo.o 625 Aso. 9375)2 with the antimony substi¬ 
tution. “x %” denotes that the length becomes large by x percent from the 
experimental value. 


Anion height (A) Mean As-Fe-As angle cr (°) 
Experimental structure^^ 1.418149 108.079 

Optimized structure 1.26538 (-10.77 %) 114.36483 (-1-5.82 %) 

Sb in chainlike As layer 1.24972 (-11.88 %) 115.380 (-^6.76 %) 

Sb in FeAs layer 1.283721 (-9.48 %) 113.9733 (+5A5 %) 



Fig. 2. (Color online) Fermi surfaces for the parent compound CaFeAs 2 . 


Ref. 15. The difference from the Fermi surfaces shown in 
Ref. 15 is the position of the three-dimensional ellipsoidal 
Fermi surface. In their calculation, there is the ellipsoidal 
Fermi surface around Z = (0,0,1 /2) point. 

Next, let us discuss antimony doping effects on the elec¬ 
tronic structure. To calculate the band structure of the par¬ 
ent compound CaFeAs 2 , we use the structural parameters ob¬ 
tained by the experiment in Ref. 5. Note that we do not use 
the optimized structure, since the lattice parameters obtained 
by the optimization are much differ from those obtained by 
the experiment. Thus, the optimization technique does not al¬ 
low to evaluate accurately the electronic structure. To assess 
the electronic structure of the antimony doped compound, we 
modify the the structural parameters of the parent compound 
obtained by the experiment. With fixing the ratios of a^lash, 
bolbsh, and /zo/^sb^ we adopt the lattice parameters of 

the 2 X 2 X 1 supercell of CaFeAs 2 are a = 7.939243 A, 
b = 7.810799 A and c = 10.317950 A. Here, cio(sh)^ ^o(Sb) 
and co(sb) are the lengths of the unit cell, and /zo(sb) is the 
anion height from Fe layer in parent (Sb-doped) compound, 
respectively. This treatment with fixing ratios of the lattice 
parameters and the anion height was used in the calculation 
to discuss the pressure effect in other pnictides.^^^ The band 
structures are obtained in folded Brillouin zone of the su- 
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Table I. Calculated crystal structure of CaFe(Sbo.o 625 Aso. 9375)2 with and without the antimony substitution, “x %” denotes that the length becomes large by 
X percent from the experimental value. 



a (A) 

b(A) 

c(A) 

a(°) 

AO 

rO 

Experimental structure^ ^ 

3.94710(10) 

3.87240(10) 

10.3210(7) 

90 

91.415(6) 

90 

Optimized structure 

3.955536 (-^0.21 %) 

3.893495 (-^0.54 %) 

10.091231 (-2.23 %) 

90.000000 

91.184438 

90.000000 

Sb in chainlike As layer 

3.978106 (-fO.786 %) 

3.9266745 (-^1.40%) 

10.088251 (-2.26 %) 

90.001299 

91.067423 

89.998870 

Sb in FeAs layer 

3.9837315 (-^0.93 %) 

3.9137245 (-^1.07 %) 

10.120249 (-1.95 %) 

90.000294 

91.257985 

90.002559 


percell. We show the unfolded and folded Brillouin zones at 

= 0 in Fig. 3. The symmetric points in the triclinic Bril¬ 
louin zone are F = (0,0,0), B = (1/2,0,0), F = (0,1/2,0) 
and G = (0,0,1/2). In this Brillouin zone, there are the very- 
narrow cylindrical Fermi surfaces on the F-F line. In the case 
of the antimony substitution with the use of the above tech¬ 
nique, the band along G-F-line shifts, where the size of the 
three-dimensional Fermi surface becomes relatively large as 
shown in Fig. 4. These results suggest that the calculated band 
structures with and without the antimony doping are similar 
to each other. In terms of the enhancement of Tc with the an¬ 
timony doping, there are three possibilities of the role of an¬ 
timony doping effects as follows. First, there is a possibility 
that these small change of the band structure brings the large 
increase of Tc. A theoretical microscopic calculation is re¬ 
quired, such as solving the Eliashberg equations. The second 
possibility is that the enhancement in this material is not 
related to the change of the band structure. One might have 
to consider the local change of the electronic states which is 
hard to be treated by the band calculation and/or the possibil¬ 
ity of the increase of the electron-electron interaction due to 
the antimony doping. The third one is the possibility that the 
adopted lattice parameters with fixing the ratios of the struc¬ 
tural parameters is not suitable to describe the electric struc¬ 
ture with and without the antimony doping. It is important to 
find out the appropriate structure which can describe the elec¬ 
tric structure with the antimony doping correctly. It remains 
as a future work. 

Finally, let us comment on the possible role of the chainlike 
As layers. Along the F-F line, there is a steep band intersect¬ 
ing the Fermi level, which originates from the As-p orbitals 
of the chainlike layers as shown in Fig. 5. This means that the 
chainlike As layers in this material are not insulating. This is 
unusual compared to other iron-based superconductors, where 
the superconducting transition temperature becomes higher 
when the “insulating” layer between the FeAs layers is more 
insulating. This also might be of some importance to under¬ 
stand the relatively high of CaFeAs 2 . 

In summary, we studied the antimony doping effect in the 
iron-based superconductor CaFe(SbjcAsi_jc )2 with the use of 
the first-principle calculation. The calculations revealed that 
the substitution of the doped antimony atom into the chainlike 
As layers is more stable than that in FeAs layers. This predic¬ 
tion can be checked by experiments. Our results suggested 
that the chainlike As layers, which exist only in novel 112 
system, have crucial role for increasing the critical tempera- 
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Eig. 3. (Color online) Schematic figure of the unfolded and folded Bril¬ 
louin zones at = 0. 
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Eig. 4. (Color online) Electronic structures with and without the antimony 
substitution onto the chainlike As layer. The experimental structural parame¬ 
ters are used for the parent compound and the modified ones are used for the 
doped compound. 


ture. With the antimony substitution, the band along G-F-line 
shifts, where the size of the three-dimensional Fermi surface 
becomes large. Since the calculated band structures with and 
without the antimony doping are similar to each other, more 
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theoretical efforts will be necessary to clarify the origin of the 
Tc increase upon antimony substitution. 
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